Neutrinos in IceCube/KM3NeT as probes of Dark Matter Substructures in Galaxy Clusters 
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Galaxy clusters are one of the most promising candidate sites for dark matter annihilation. We 
focus on dark matter (x) with mass in the range (lOGeV — 100 TeV), annihilating through the 
channels \X ~* XX XX ~^ tt, or \X ~ * vvvv, and forecast the expected sensitivity to 

the annihilation cross section into these channels by observing galaxy clusters at IceCube/KM3NeT. 
Optimistically, the presence of dark matter substructures in galaxy clusters is predicted to enhance 
the signal by (2 — 3) orders of magnitude over the contribution from the smooth component of 
the dark matter distribution. Optimizing for the angular size of the region of interest for galaxy 
clusters, the sensitivity to the annihilation cross section, (cry), of heavy DM with mass in the range 
(300GeV - lOOTeV) will be O(10" 24 cn^s" 1 ), for full IceCube/KM3NeT live time of 10 years, 
which is about one order of magnitude better than the best limit that can be obtained by observing 
the Milky Way halo. We find that neutrinos from cosmic ray interactions in the galaxy cluster, 
in addition to the atmospheric neutrinos, are a source of background. We show that significant 
improvement in the experimental sensitivity can be achieved for lower DM masses in the range 
(10 GeV — 300 GeV) if neutrino-induced cascades can be reconstructed to « 5° accuracy, as may be 
possible in KM3NeT. We therefore propose that a low-energy extension "KM3NeT-Core" , similar 
to DeepCore in IceCube, be considered for an extended reach at low DM masses. 
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I. INTRODUCTION 

There is overwhelming evidence for, yet unexplained, 
invisible mass in our Universe [1-4]. Particles in the 
standard model of particle physics cannot account for the 
major fraction of this excess mass, but a new particle 
with weak-scale annihilation cross sections to standard 
model particles, as predicted in several extensions of 
the standard model of particle physics, would naturally 
explain its observed abundance [5-9] . This has motivated 
a comprehensive search for the particle identity of this 
"dark matter" (DM) using (i) direct production of DM 
at colliders [10, 11], (ii) direct detection of DM via 
elastic scattering [12-23] and (hi) indirect detection of 
DM via its annihilation or decay [24-44]. This three- 
pronged approach to DM detection is necessary because 
a single experiment cannot probe all the properties of 
DM. For example, collider experiments mainly probe 
production of DM particles, whereas direct detection only 
probes the interaction between the DM particle and the 
particular detector material [45]. Analogously, in an 
indirect detection experiment, we learn about the final 
states of DM annihilation or decay. 

Indirect detection experiments are also sensitive to the 
DM density distribution at cosmological scales in this 
Universe unlike direct detection experiments which are 
only sensitive to the local DM distribution in the Milky 
Way. If LHC detects a DM candidate, then indirect 
detection experiments are also useful to determine 
whether that particular DM candidate makes up most 
of the DM in the Universe [46]. Indirect detection 
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experiments, looking for products of DM annihilation 
in astrophysical sources, only detect a handful of the 
final states, e.g., photons, electrons, protons, neutrinos, 
and their antiparticles. If these experiments detect 
a signal that requires a cross section larger than the 
thermal relic annihilation cross section it would challenge 
a simple thermal WIMP paradigm of DM, and thus 
provide a crucial test of the WIMP paradigm [47]. On 
the other hand, there is no guarantee that a signal must 
be found if we can probe cross sections comparable to, or 
smaller than, the thermal relic annihilation cross section 
- annihilations could proceed to undetected channels. 
In that case, however, one sets an upper bound on the 
partial annihilation cross sections into these observed 
channels, constraining particle physics models of DM. 

Several astrophysical targets, e.g., the Sun, the Milky 
Way, dwarf galaxies, and galaxy clusters, may be 
observed by indirect detection experiments. A careful 
estimate of the signal and the background for each of 
these source classes is needed to determine which of 
these targets provides the best signal-to-noise ratio for 
a given DM model. The Sun accumulates DM particles 
while moving through the DM halo of the Milky Way. 
Due to the high density at the core of the Sun, for DM 
mass > 300 GeV, annihilations products are absorbed 
and the sensitivity of DM annihilation searches weakens 
considerably, making it inefficient for probing high mass 
DM [48-52]. The Milky Way is dominated by DM in its 
central regions, but unknown astrophysical backgrounds 
make it difficult to disentangle the signal [53-58] , whereas 
the diffuse component of the Milky Way DM halo [33, 37] 
leads to a significantly reduced signal. Dwarf galaxies 
have a high mass-to-light ratio and are one of the ideal 
targets for detecting DM in gamma-ray experiments with 
subdegree angular resolution [27, 28, 59-66]. 
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Galaxy clusters have the largest amount of DM 
amongst all known classes of gravitationally bound 
objects in the Universe. Although the background 
due to other astrophysical sources is also large therein, 
the contribution of DM substructures can enhance the 
DM annihilation signal from the smooth component, 
typically modeled using a Navarro, Frenk, and White 
(NFW) profile [67]. This enhancement depends 
on the abundance of DM substructures. State-of- 
the-art galaxy cluster simulations do not have the 
resolution to directly calculate the contribution due to 
the theoretically expected least massive substructures. 
However using theoretically well-motivated values for 
the mass of the smallest substructure and extrapolating 
the abundance of substructures to these lowest masses, 
high resolution computer simulations predict that galaxy 
clusters provide the best signal-to-noise ratio for DM 
annihilation signal [68]. Note that even a moderate 
enhancement due to DM substructure, as advocated 
in [65] following the works in [69, 70] predicts that galaxy 
clusters give the best signal-to-noise ratio for analysis 
where the field-of-view is greater than or equal to 1°. 
This strongly motivates observations of galaxy clusters 
to search for DM annihilation signals [71-79]. 

Neutrino searches, among other indirect searches for 
DM, have distinct advantages. Being electrically neutral 
and weakly interacting, neutrinos travel undeflected 
and unattenuated from their sources. So neutrinos 
can provide information about dense sources, which 
may be at cosmological distances, from which no other 
standard model particles can reach us. Another crucial 
motivation to look for neutrinos is that many standard 
model particles eventually decay to produce neutrinos 
and gamma rays as final states. Detecting neutrinos is 
therefore complementary to gamma ray searches from 
DM annihilation, which have become very exciting in 
recent times [27, 29, 30, 33, 35]. For very heavy DM, the 
gamma rays produced in the DM annihilations cascade 
and the constraints on DM annihilation cross sections 
become weaker than those obtained using neutrinos. 
Also, for hadronic explanations of any gamma ray 
and cosmic ray excesses, detecting neutrinos will be 
a smoking gun signature. Finally, direct annihilation 
to neutrinos is impossible to detect using any other 
detection channel, with electroweak bremsstrahlung 
being a notable exception [80] although the limits 
obtained in that case turns out to be weaker than those 
obtained by direct observation of neutrinos [81]. In fact, 
neutrinos, being least detectable, define a conservative 
upper bound on the DM annihilation cross section to 
standard model particles [82, 83]. 

Limits obtained by gamma ray telescopes are typically 
stronger than that obtained using neutrino telescopes, 
but the larger angular resolution of a neutrino telescope, 
compared to a gamma ray telescope, means that the 
results obtained in a neutrino telescope is less dependent 
on the central part of the DM density profile (which 
gives the strongest signal in a gamma ray telescope) 



where the uncertainty obtained in DM simulations is the 
largest. Neutrinos telescopes are also able to view a 
target source for a longer time compared to a gamma 
ray telescope, though this advantage is mitigated by 
the smaller cross section of neutrino detection. Another 
advantage of neutrino telescopes is that they are able to 
view a large number of sources simultaneously and can 
be used to find dark matter in a region which is dark in 
the electromagnetic spectrum. These arguments and the 
availability of large neutrino telescopes strongly motivate 
a search for DM annihilation using neutrinos. 

Although dwarf galaxies are known to be the 
best targets for dark matter searches for gamma-ray 
experiments, they are not the best targets for neutrino 
experiments. The reason for this is the limited angular 
resolution of a neutrino telescope, which is > 1°. Dwarf 
galaxies have an angular size of < 1° and thus when 
a neutrino telescope takes data from a dwarf galaxy, 
even with the minimum angular resolution, the size of 
the dwarf galaxy is smaller than the data-taking region, 
which implies a worse signal-to-noise ratio. However, 
galaxy clusters have a typical size of a few degrees and 
hence even when neutrino telescopes are taking data in 
the larger than minimum angular resolution mode, the 
size of the galaxy cluster fills up the entire data-taking 
region. This ensures that, unlike in the case of dwarf 
galaxies, there is no position in the data-taking region 
from where there is no potential signal candidate and 
thus provides a better signal-to-noise ratio. 

Neutrinos from galaxy clusters have been considered 
previously by Yuan et al. [84]. In that paper, the DM 
halo for a galaxy cluster was obtained from extrapolation 
of the DM halo obtained from the simulation of a 
Milky Way like galaxy [85, 86]. Using the Fcrmi- 
LAT limits from galaxy clusters, Yuan et al. constrained 
the minimum DM substructure mass, and analyzed 
muon tracks in IceCube to obtain a constraint on DM 
annihilation cross section. 

In this paper, we investigate neutrinos from galaxy 
clusters using the latest DM density profiles, as given in 
Gao et al. [68]. This gives us updated inputs for both 
the smooth and the substructure components of DM 
in galaxy clusters. For comparison, we also calculate 
our results by taking the smooth and the substructure 
components of DM profile from the work by Sanchez- 
Conde et al. [65] and find that due to the smaller boost 
factors (about a factor of 20 smaller boost factors than 
compared to that in [68] ) , the sensitivity of the neutrino 
telescope for this parametrization of the DM profile 
is about a factor of 20 worse than what is obtained 
while using the DM substructure modeling of [68]. 
We also take into account neutrinos produced due to 
cosmic ray interactions in the galaxy cluster, ignored in 
previous studies. With these updated inputs, we analyze 
the expected signals and backgrounds at IceCube and 
KM3NeT for both track and cascade events. While, 
quantitative improvement in the detection prospects is 
found for track searches, qualitative improvement in 



sensitivity and reach at low DM masses is expected if 
KM3NeT deploys a low energy extension, which we call 
KM3NeT-Core, and is able to reconstruct cascades with 
a pointing accuracy down to 5° as claimed by Auer [87] . 

The remainder of the paper is arranged as follows. 
In Sec. II, we discuss the neutrino flux from DM 
annihilation, using the DM density profile of a typical 
galaxy cluster. In Sec. Ill, we discuss neutrino detection 
and relevant backgrounds at a neutrino telescope. In 
Sec. IV we discuss the results, showing our forecasted 
sensitivity to (av) for the considered annihilation 
channels, and conclude in Sec. V. 



II. DM DISTRIBUTION AND NEUTRINO 
PRODUCTION IN GALAXY CLUSTERS 

The number flux of neutrinos per unit energy interval 
(in GeV _1 cm _2 s _1 ) for a given final state of DM 
annihilation is given by [78] 

where m x denotes the mass of the DM particle (in units 
of GeV), (av) denotes the thermal-averaged annihilation 
cross section into the final state which can produce 
neutrinos (in units of cm 3 s _1 ). dN u /dE l , denotes the 
energy spectrum of the neutrinos from the various final 
states of DM annihilation (in units of GeV -1 ). The 
integral J dl p 2 [r(l),ip] is the line-of-sight integral of the 
DM density distribution, with / denoting the line-of-sight 
distance (in units of cm), p(r) denoting the DM density 
distribution function at a point r (in units of GeVcm -3 ). 

We have assumed here that the DM is its own 
anti-particle, which gives an extra factor of 2 in the 
denominator of the expression in Eq. (I). We also 
assume that the galaxy cluster is close enough so that 
the neutrino energy is not red-shifted significantly. We 
emphasize that even for a neutrino telescope, a nearby 
galaxy cluster is not a point source and hence an angular 
dependence, tp, of the line-of-sight integral is present. 
Therefore, we have to integrate over the relevant solid 
angle, Ail = 2tt J ma * simp dip, where f/'max is the 
angular radius of the region of interest. 

It can be seen that the neutrino flux is written in such 
a way that it is a product of the astrophysics quantities, 
J dl p 2 [r(l),tp], with the particle physics quantities, 
to~ 2 ((av) /2)dN v I dE v . In the following subsections, we 
outline how we have calculated each of these quantities 
for our analysis. 



A. DM distribution 

In this section, we describe the DM density distribution 
in a typical galaxy cluster. Although we shall refer to the 
Virgo galaxy cluster for specific quantitative details, the 
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FIG. 1. Figure to illustrate the line-of-sight integral. The blue 
shaded region is the galaxy cluster with C as its centre. The 
position of the observer is marked by the point O. The virial 
radius of the galaxy cluster is denoted by r v i r . The distance 
of the observer, O, to the centre of the galaxy cluster, C, 
is denoted by D. The distance of a typical point inside the 
galaxy cluster from the centre of the galaxy cluster and the 
observer is denoted by r and I respectively. 



same physical description is qualitatively applicable to 
other galaxy clusters. 

Galaxy clusters are the most massive gravitationally 
bound objects in the Universe today [89, 90]. A typical 
galaxy cluster has a mass of ~ 0(1O 14 M0) and is 
virialized up to a radius of ~ O(Mpc). We take the 
smooth component of the DM density profile to be 
parametrized by an NFW profile [67] 




where r s is the scale radius and p s is the DM density at 
distance ~ 0(r 8 ) from the centre of the galaxy cluster. 

Given the redshift, z, and the virial mass, M v ; r , of a 
galaxy cluster, the virial radius, be determined 

from the following relationship, as given by Ando and 
Nagai [78], 

M vir = ^7rr 3 ir A vir (z)p c (z) . (3) 

Here, virial quantities are identified by using "vir" in 
the subscript, p c (z) is the critical density of the Universe 
and the cosmological factor A v j r (z) = 82d — 39<i 2 + 18-7T 2 , 

where d = — Q\/ (fi A + fi m (1 + z) 3 ) [91]. We assume a 
ACDM model for the Universe for all calculations: fi A = 
0.73, Q m = 0.27 and H = 73 km s -1 Mpc" 1 [92]. The 
scale radius is obtained from the equation r s = r vir /c vir , 
where c v ; r denotes the concentration parameter which is 
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given by [78], 



above, we find that 



7.85 



(1 + z) - 71 V2 x lO^/i-iMe 



-0.081 



(4) 



To obtain p s , we equate the virial mass of the cluster, 
M vir , to the volume integral of p(r) up to r vir . 
Analytically, we obtain 



C virPc 



3 log(l + c vir ) 



(5) 



1 + c v 



For the Virgo galaxy cluster, the virial mass is taken 
to be M vir = 6.9 x 10 14 M & [93] and the redshift is 
taken to be z = 0.0036 [94]. Note that the redshift is 
too small to affect neutrino energies appreciably. Using 
the value of the critical density of the Universe, p c = 
0.54 x 10~ 5 GeV cm~ 3 [92], we get the virial radius of the 
Virgo galaxy cluster to be r vir — 2.29 Mpc. We use the 
concentration parameter, c vlr — 4.98, as in [78], which 
gives r s — 0.46 Mpc. For the Virgo galaxy cluster, we find 
that p s = 2.19 x 10~ 2 GcVcm~ 3 . We note that this value 
of the central DM density is about a factor of 2 lower 
than what we would have obtained if we had followed 
the prescription in Han et al. [88] . This difference can be 
traced to the fact that the cosmological factor, A v i r (2;), in 
our calculation has a value of ~ 100, whereas the similar 
expression for p s , as given in [88], gives the cosmological 
factor to be 200. We adopt the optimistic value of A v j r , 
and hence use p s — 4.38 x 10~ 2 GeV cm -3 throughout 
this work, but remind the readers that a lower value of 
A v ; r (,z) by a factor of two can decrease the annihilation 
signal and sensitivity by a factor of 4. 

In Fig. 1, we schematically show how to calculate the 
line-of-sight integral. Here O is the position of the 
observer and the blue shaded region is the galaxy cluster 
whose centre is denoted by C. The virial radius of the 
galaxy cluster is shown as r V j r and the distance to the 
centre of the galaxy cluster is denoted by D. The line-of- 
sight distance to a point inside the galaxy cluster which 
is at a distance r from the centre of the galaxy cluster is 
given by I. 

The line-of-sight integral, as a function of the angle ip, 
is defined as 



j(V) = j i dlp 2 {r(l),^], 



(6) 



where 

r = y/P + D 2 - 2DI cos ip , 
Z maJ c,min = £>cosV ± J D 2 cos 2 ijj - (D 2 - rV) . (8) 



(7) 



We call this integral the j-factor for future reference. 
The distance to Virgo galaxy cluster is taken to be 
D = 19.4 Mpc [94]. Using the parameters mentioned 
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dnj NFVf (ip) 



= 2.064 x 10~ b GeV 2 cm~ b Mpc , (9) 

where Afi = 2tt Jq ax simp dip, and i/'max w 6° for the 
Virgo galaxy cluster. Recent high resolution simulations 
of galaxy clusters, in particular the Phoenix project [68], 
show a high concentration of DM substructures in 
addition to the smooth NFW profile. Tidal forces destroy 
the smallest mass substructures in the inner regions of 
the galaxy cluster so the inner region of a galaxy cluster 
(< 1 kpc) is dominated by the smooth NFW profile. 
However, the DM density in the outer region of a galaxy 
cluster is dominated by the DM substructures [68] . This 
suggests that one should search for extended emission 
while looking for DM annihilation signal from a galaxy 
cluster. 

The contribution to the DM annihilation due to 
substructures depends on their abundance. Recent 
simulations can only resolve substructures of masses 
> 10 5 M© but theoretical considerations suggest that 
the minimum substructure mass for cold DM is in the 
range (10 -12 M Q - 1O" 6 M ) [95]. In order to obtain 
the DM annihilation signal, we have to extrapolate the 
substructure abundance, using a halo mass distribution 
function from the simulations, from a mass of ~ 10 5 M Q 
to a minimum substructure mass of ~ lO~ 6 A/0. This 
11 orders of magnitude extrapolation is the largest 
source of uncertainty in our calculation. However, it 
must be emphasized that even with a mass resolution 
of ~ 5 x 10 7 M Q , the galaxy cluster simulations predict 
that the substructure contribution completely dominates 
the smooth contribution at radii > 400 kpc [68]. 

Assuming the smallest substructures to have masses 
~ 10~ 6 M Q , Han et al. [88] parametrize the j-factor due 
to substructures as 



jsub(^) 



b(M 2 oo)Ji 
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Jsub(^) 



jsubCV^oo) e 



-2.377 



(V74) 2 



, ^200 , 



(10) 



(11) 



where 6(M 200 ) = 1.6 x 10- 3 (M 200 /M Q ) - 39 is the boost 
factor. Here M200, ^200, and ?"2oo are the mass, angular 
radius, and radius of the cluster where the average DM 
density is 200 times the critical density of the Universe. 
Using the value of M 2 oo, as given in [88], we obtain the 
boost factor, &(M 2 oo) ~980. As mentioned in [68], this 
boost factor is about an order of magnitude larger than 
the analogous boost factor obtained from galaxy halos. 
A boost factor of ~ 1000 for galaxy clusters was also 
analytically obtained in [76]. 

Here we again mention that if we follow the galaxy 
cluster DM substructure modeling of [65], the boost 
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factor that we obtain is 55 for the Virgo galaxy cluster 
and between 34 and 54 for other galaxy clusters that 
were considered in that work. Hence there is a factor 
of ^20 uncertainty in the sensitivity to DM particle 
properties that can be derived from observation of galaxy 
clusters both by gamma-rays observations and neutrino 
observations. 

We scale the line-of-sight integral j(ifj) to our local DM 
density-squared and to the distance to the Galactic centre 
from the Sun to define the scaled j-factor 



dl ( P [r{l)M 
.5kpc V0.3GeVcm- 3 



(12) 



In Fig. 2, we plot against angle if), for the Virgo 

galaxy cluster for both the DM profile models in [65, 88]. 
It is easily seen that the presence of substructure provides 
a large boost to the DM annihilation signal for both 
the DM profile models, although the boost factors are 
different for both the models. The contribution from 
the NFW halo is concentrated at the centre whereas the 
contribution from the DM substructure is more extended 
for both the DM profile models. We use the model in [88] 
for all our subsequent results. To obtain the results 
for the DM profile modeling with [65], one can simply 
decrease the sensitivity in the result section IV by a factor 
of ^20. We remind the reader that at present due to 
the limited numerical resolution of the DM simulations, 
it is impossible to completely resolve the question of 
the boost factor which not only depends on the lowest 
DM substructure mass but also on merging of different 
galaxies to form a galaxy cluster. 
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FIG. 2. Scaled line-of-sight integral (l.o.s.), as a 

function of the angle ip (see Eq. 12). The black line shows 
the scaled l.o.s. for an NFW halo. The red dot-dashed line 
represents the scaled l.o.s. due to DM substructures following 
the work of [88]. We also show the scaled l.o.s. due to 
DM substructures following the work of [65] by the brown 
dash-dash-gap line. The blue dotted line is the combined 
contribution of the NFW halo and the substructure following 
the work of [88] . The combined contribution of the NFW halo 
and the substructure following the work of [65] is also shown 
by the green dash-dot-dot line. H in parenthesis denotes 
parametrization taken from [88] where SC in parenthesis 
denotes parametrization taken from [65]. Figure adapted 
from [881. 



B. Neutrino spectrum at source 

Now we turn our attention to the particle physics 
relevant for calculating the neutrino flux from DM 
annihilation. Since the DM in galaxy clusters is non- 
relativistic, with typical velocities v ~ 10 3 kms _1 , the 
DM annihilation products in a 2-body final state with 
identical particles are produced with an energy equal to 
the mass of the DM particle. 

In this paper we study the the sensitivity of neutrino 
telescopes to (crv) for DM annihilation to four interesting 
channels: (i) \X M + A*~ , (u) XX ~> vV , ( ni ) XX it 
and (iv) \X ~> W ~> vvvv. All these chosen final state 
particles have or produce neutrinos on decay, and we 
forecast the sensitivity to the annihilation cross section 
that can be obtained using a neutrino telescope. 
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Neutrino oscillations ensure that there is a 1:1:1 ratio of 
the fluxes of the v e , v„ and v T reaching the detector. An 
analogous equation holds true for antineutrinos. 



XX 



1- XX "> A 4 A* 

The XX ~ * channel leads to signals in both 

gamma rays and neutrinos, and therefore quite promising 
for multi messenger studies. The normalised neutrino 
spectrum from decays of energetic muons of energy is 



Searching for direct annihilation to neutrinos is 
motivated by the presence of sharp spectral feature in 
the neutrino spectrum in the channel. Although this 
channel is suppressed for a Majorana or a scalar DM 
particle, there exist models in which the DM coupling 
to neutrinos is enhanced. This channel also gives the 
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most stringent limits for DM annihilation in a neutrino 
telescope. 

The neutrino spectrum due to direct annihilation to 
neutrinos is given by 



dN, 
dE, 



- = 5{E U - m x ) 



(15) 



Due to the finite energy resolution of the neutrino 
telescope, the dirac-delta function gets smeared out. We 
model the neutrino spectrum as a gaussian with centre at 
m x and a full- width at half-maximum given by the energy 
resolution of the detector [97]. Neutrino oscillations 
ensure that there is a 1:1:1 ratio of the fluxes of the v e , 
and v T reaching the detector. 



3. 



XX 



tt 



The third channel which we consider is \X ~* ti- This 
is the most favored annihilation channel, from helicity 
arguments, for heavy (> 175 GeV) DM, if the DM is a 
Majorana fcrmion or a scalar. 

The top quark decays to W^-boson and a 6-quark with 
a branching ratio of > 99% and the subsequent decay 
of VF-boson and hadronization of the &-quark produces 
neutrinos. As an approximation, we consider only the 
prompt neutrinos produced by the decay of the W^-boson 
and the 6-quark. Following [98], we derive the highest 
energy muon neutrino flux due to the top quark decay as 
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where I denotes the relevant lepton states in the decay 
of the T^-boson and the decay of the &-hadrons. The 
corresponding branching ratio for the decay of the W- 
boson and the 6-hadrons is denoted by T, and the 
corresponding values are taken from PDG [92]. The 
Lorentz factor is denoted by 7t = E t /m t = l/y/l — {3?- 
Ew and fiw are the energy and velocity of the VK-boson 
in the top quark rest frame. E± = E v ^ x /(l =p (3 t ) 
represents the maximum and minimum energy of the 
neutrino in the moving frame of the top quark. The 
limits of the neutrino energy in the moving frame of 
the W-boson is denoted by e± = E w (l ± p w )/2. If 
the energy of the 6-quark in the rest frame of the top 
quark is denoted by £7, then the hadronization energy 
can be approximated as Ed = ZfEb where we take 
the value of zj from [98]. The function Db[x,y] — 
I (6 ln(y/x) + 4(y 3 - x 3 ) + 9(x 2 - y 2 )). We ignore the 
lower energy muon neutrinos produced due to the decay 
of the muons in the final state. Neutrino oscillations 



ensure that there is a 1:1:1 ratio of the fluxes of the v ei 
and v r reaching the detector. 



4- XX 



VV 



This channel is motivated by the secluded DM 
models [99, 100], in which the DM annihilates to two 
light vector bosons V (or a similar mediator) each of 
which then decay to standard model particles and can 
be observed [50, 101, 102]. If the decay is primarily to 
neutrinos, one gets two neutrino pairs in the final state. 
There is a recent proposal in [103], which addresses some 
of the purported small-scale problems in ACDM, also the 
DM annihilation to neutrinos is enhanced and hence this 
model can be tested using neutrino telescopes. 

The neutrino spectrum has a box-like structure 



dN u _ 4 
~d~E~ v ~ ~AE 



Q(E U - E-)Q(E+ - E v 



(17) 



where denotes the Heaviside-theta function. The 
maximum and minimum energy of the neutrino in this 

case is denoted by E± = (m x ± Jm 2 , — rriy)/2. The 



width of the box function by AE = W m x — m\ 

Neutrino oscillation ensures that the ratio of the neutrino 
fluxes reaching the neutrino detector is 1:1:1. 



III. DETECTION AND BACKGROUNDS 



Neutrino detection 



In a km 3 -scale neutrino telescope like IceCube [104] or 
KM3NeT [105], neutrinos are detected as two different 
types of events: tracks and cascades. 



1. Tracks 

The tracks are produced by the charged current 
interaction of the muon neutrinos and antineutrinos. At 
these high energies, the muons are produced by muon 
neutrinos interacting with the detector material or with 
the surrounding medium and the muon track is generally 
not contained inside the detector [106]. Due to the long 
range of the muon tracks, the effective volume of the 
detector is increased and the increase in the volume is 
determined by the range of the muon, of a given energy 
E, given by integrating the energy loss rate 



dE 
dX 



= a + /3E , 



(18) 



where X denotes the column density (in units of gem 2 ). 
For our calculations, we take a = 2MeVcm 2 g _1 and 
= 4.2 x 10- 6 cm 2 g- x [106]. 
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The number of neutrinos detected per unit energy 
interval for muon tracks, which are not contained inside 
the detector is given by [106] 



dE„ 



tracks 



N AP T A det 



dE„ : ^<Jcc{Ev)e 
dE v 



(19) 



In the above formula, N A denotes the Avogadro's 
number, p represents the density of the detector material, 
T is the time of observation, ^4dot denotes the area of 
the detector, acc{E v ) denotes the charged current cross 
section of the muon neutrino with the detector material 
or its surroundings, L the length traveled by the neutrino 
in the Earth, and A is the mean free path of the neutrino. 

The factor A^p^ 1 / (a + f3E^) accounts for the 
increased volume of the detector due to the long muon 
range. We take A^ct = 1 km 2 and T = 10 years for the 
XX -> XX -> vv, XX -> it, and \X ~> vvvv 

channels. The values of <Jcc(E v ) are taken from [107]. 
The exponential suppression is due to the absorption 
of very high energy neutrinos (> 100 TeV) as it passes 
through the Earth. The mean free path of the neutrinos 
in Earth matter is given by A = l/(n<7 to t), where n 
denotes the number density of target particles and a to t 
denotes the total neutrino-nucleon cross section, which 
we take from [107]. For the energies considered here, 
the exponential factor is ~ 1. 

For non-contained muon track events in IceCube, the 
energy is obtained by using Eq. (18) after measuring the 
muon energy loss inside the detector [108]. The limits of 
the integral in Eq. (19) imply that a muon of energy E^ 
can be produced by any with an energy > E^. 

The energy range that we explore using muon tracks 
is (100 GeV - 100 TeV). Energy resolution for muon 
tracks is approximately 0.3 in units of logio-E [97]. 
Following [109], we take the energy bin for signal 
calculation to be (max(_E t j ires , m x /5), m x ). This energy 
bin is much bigger than the energy resolution of 
IceCube/KM3NeT [108]. We expect a full spectral 
analysis by IceCube/KM3NeT collaboration to give 
much better sensitivity as the shape of the signal and 
background spectra are very different. In this regard, 
the results presented here can be treated as conservative. 

Angular pointing for tracks is quite accurate. For 
neutrino energies > 100 GeV, the angular resolution is 
within 0.5° and 1° [108]. 



to the charged current interaction of the v T with the 
detector material produces a hadronic cascade from its 
hadronic decay products and an electromagnetic cascade 
from the electrons arising from r decay. The non- 
leptonic final states in a neutral current interaction 
causes a hadronic cascade in the neutrino telescope. 
These cascades are contained inside the detector, act 
as almost point sources of light, and are calorimetric. 
The cascade search also has lower atmospheric neutrino 
background [110]. Cascades has been detected in 
IceCube [111] and recently also in DeepCore [112]. 

The number of neutrino events detected via cascades 
per unit energy interval is given by [109] 



dN v 



dE v 



= N A T V casc 



dE„ 



dE„ 



,(20) 



where V casc = 0.02 km 3 is the volume available for 
cascades in a detector like IceCube-DeepCore and <7nc 
the neutral current cross section of neutrinos, which we 
take from Ref. [107]. Other symbols have meanings and 
values as previously defined. 

The mass range of DM that we explore in the cascade 
analysis is (10 GeV - 10 TeV). The energy resolution 
for cascade like events is approximately 0.18 in units of 
logio-E [97]. Following [109] we take the energy bin for 
signal calculation to be (max(£ , t hres, w x /5), m x ). This 
energy bin is much larger than the energy resolution 
of IceCube/KM3NeT [108]. We expect a full spectral 
analysis by IceCube/KM3NeT collaboration to give 
much better sensitivity as the shape of the signal and 
background spectra are very different. In this regard the 
results presented here can be treated as conservative. 

Achieved angular resolution of cascades in IceCube is 
about 50°, but is expected to be significantly improved 
in the future with more advanced reconstruction 
algorithms in DeepCore [113]. With a large angular 
resolution the background due to atmospheric neutrinos 
is overwhelming, and improving the resolution drastically 
cuts down background. Encouragingly, Auer [87] 
discusses a procedure which can be used to reconstruct 
the angular resolution of cascades to about 5° in 
KM3NeT. We shall show that with such improved 
angular resolution, the sensitivity to DM annihilation 
cross section by cascades increases significantly. 



2. Cascades 

Charged current interactions of v e and v T and their 
antiparticles, and neutral current interactions of all 
flavors of neutrinos produce cascades. The electron 
produced due to the charged current interaction of the 
v e with the detector material causes an electromagnetic 
cascade in the detector. The r-lepton produced due 



B. Detector Configurations and Backgrounds 

While calculating the sensitivity to the DM 
annihilation cross section using muon tracks, we 
assume that the neutrino telescope only looks at upgoing 
tracks. This means that IceCube will look at galaxy 
clusters in the northern hemisphere and KM3NeT, 
while using muon tracks for their analysis, will look at 
galaxy clusters in the southern hemisphere. Looking 



for upgoing tracks eliminates the background caused by 
downgoing atmospheric muons. 

For the cascade analysis we shall assume that KM3NeT 
includes a DeepCore-like low energy extension, which 
we call KM3NeT-Core. We assume the mass of the 
KM3NeT-Core to be the same as that of DcepCore. 
Similar to DeepCore, we shall assume that KM3NeT- 
Core will use the remainder of the KM3NeT as a veto. 
Such an arrangement allows the low energy extension in 
KM3NeT to have a Air field of view, and therefore this 
low energy extension in KM3NeT can also be used to 
detect galaxy clusters in the northern hemisphere. 

With the configurations explained in the previous two 
paragraphs, the backgrounds in both track and cascade 
analyses are due to atmospheric neutrinos and neutrinos 
from cosmic ray interactions in the galaxy cluster. 
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The measured atmospheric za. 



flux for E,. in the 



range of (100 GeV - 400 TeV) is reported in [108]. 
The measured spectrum is fit well by the angle-averaged 
atmospheric neutrino spectrum given in [114]: 



d<j> a 



atm ^„ 1 x 10" GeV^km- V -1 (21) 
In (1 + 0.024 £„) In (1 + 0.00139 



1.33 



0.201 E v 



where 3>atm = 1*95 f° r neutrinos and 1.35 for 
antineutrinos, and E v is in GeV. The atmospheric v e 
flux is taken from [115] and the v T flux is from [116]. 

In addition to the atmospheric neutrino, the neutrinos 
produced by cosmic ray interactions inside the galaxy 
cluster also acts as an additional background. We take 
the neutrino flux produced in cosmic ray interaction 
in the galaxy clusters from the calculations by 
Murase et al. [117]. They consider acceleration of cosmic 
rays with energies between 10 17 eV and 10 185 eV in 
shocks in galaxy clusters. In a 1° x 1° angular bin, 
they estimate <1 (u^ + V^) event per year above 1 TeV. 
Although this estimate is somewhat model-dependent, 
we emphasize that this is an essential background that 
one has to take into account while searching for neutrinos 
from DM annihilation in galaxy clusters. If it turns out 
that galaxy clusters are not sources of cosmic rays in 
this energy range then this background can be lower but 
we assume the larger background rate for conservative 
estimates. 



IV. RESULTS 

In this section we calculate the neutrino fluxes 
observed for the four annihilation channels chosen above, 
and compare them with the expected backgrounds to 
determine the sensitivity in the (av)-m x plane for each 
channel. However, before we proceed to results specific 
to each channel, we identify some broad features. 

The first key result is regarding the optimal size of 
the region of interest. The signal we are looking for 
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FIG. 3. Ratio of approximate signal-to-noise versus the 
angular size of the chosen region of interest around the Virgo 
galaxy cluster. The vertical lines show the values of some 
representative angular radii in degrees. 



is proportional to (j su b + Jnfw) dtt. We scale this 
quantity with the local DM density squared and the 
distance to the Galactic centre from the Sun, as in 
Eq. (12), to get 

/■An 

^tot(V'max) = / (jsub + J*NFw) dCl , (22) 

Jo 

where Af2 depends on the angular radius ipmax of the 
region of interest. The fluctuations in the atmospheric 
neutrino background are proportional to yj A57(?/> max ). 
Therefore, the signal-to-noise ratio is approximately 
proportional to Jtot(V'max) / V Ail^ max ) ■ 

In Fig. 3 we show this ratio for the Virgo galaxy cluster, 
as a function of the chosen angular radius V'max of the 
region of interest. We can see that due to the extended 
nature of the DM substructure profile, a region with 
s» 2° angular radius around the galaxy cluster gives the 
best signal-to-noise ratio. We have verified using our 
numerical results that the sensitivity obtained with a 2° 
observation window is about a factor of 1.5 better using 
a 1° window. 

Neutrino telescopes should therefore carefully optimize 
for the observation window. Selecting a circular region of 
diameter ~ 4° around the centre of the galaxy cluster and 
accepting signal events coming from that circular region 
appears to give the best signal-to-noise ratio. Depending 
on the specific DM profile of a galaxy cluster, this choice 
of angle may change but, in general, we expect that, for 
any nearby galaxy cluster the best signal-to-noise ratio is 
achieved by treating it as an extended source, as opposed 
to a point source. 
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FIG. 4. Sensitivity to (av) versus the DM mass, for the annihilation channels XX /- t+ A t_ (t°P panels), \X ~^ v ^ (second panel 
from top), XX ~> ti (third panel from top), and \X ~ * vvvv (bottom panel) using 10 years of IceCube/KM3NeT-Core data for 
all the channels. Left panels: Sensitivity that can be obtained from muon tracks. Right panels: Sensitivity that can be obtained 
using cascades. The lines are labeled with the assumed DM density profile: smooth component only (NFW) or including both 
smooth and substructure components (NFW+sub), as well as the angular radius chosen for the region of interest (1°, 2°, 5°, 
and 10°). We use the DM substructure modeling of [88]. A more conservative DM substructure modeling following [65] will 
worsen the sensitivity in each case by a factor of ~20. 
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The second key result is that the presence of 
substructures gives three orders of magnitude more 
promising results than the smooth NFW profile alone. 
This boost provided by the substructures make galaxy 
clusters an exciting target for neutrino telescopes. Using 
track-like events, the sensitivity is typically in the range 
(av) > (10~ 24 -10~ 22 ) cm 3 s" 1 at DM masses (100 GeV- 
100 TeV). At lower masses, the sensitivity worsens 
quickly because the events are below threshold (note 
the upturn in Figs. 4 and 5). However the sensitivity 
also worsens at extremely high masses because as m x 
increases, the number of DM particles decrease for a given 
DM density, reducing the annihilation fluxes. 

The third key result we find is that if KM3NeT can 
reconstruct cascades with an angular resolution of ~ 5° 
and has a DeepCore-like low energy extension, a new 
window of observation opens up at lower DM masses 
(10 GeV — 100 GeV). The sensitivity of neutrino cascade 
observations remain competitive with track analyses at 
masses up to 10 TeV. This complementary measurement 
of muon tracks and cascades may be useful to explicitly 
determine the neutrino flavors in the DM annihilation 
products. We believe this should encourage the KM3NcT 
collaboration to improve their cascade pointing to < 5° 
and include a DeepCore-like low energy extension in 
KM3NeT. 

We emphasize that if neutrino telescopes detect a DM 
annihilation signal from galaxy clusters at a sensitivity 
forecasted in this paper and if the emission profile is 
found to be extended, it will be a strong indication for 
the presence of substructures. If the neutrino signal 
is not extended but the cross section is comparable 
to what is forecasted to be testable then it will favor 
an enhanced annihilation cross section, rather than the 
presence of substructures. A particle physics explanation 
of the enhanced DM annihilation cross section will 
then be required [100, 118, 119]. If a signal is not 
detected at an annihilation cross section testable at 
neutrino telescopes, then it will either constrain the 
minimum DM substructure mass and abundance or the 
annihilation cross section. In that situation, one will 
have to first infer the DM annihilation cross section 
from some other astrophysical source to infer something 
about the minimum DM substructure mass and the DM 
substructure distribution. 



A. XX -»■ MM 

In Fig. 4 (top left panel), we show the sensitivity to 
(av) for the DM annihilation to M + M~i that can be 
achieved by observing muon tracks at IceCube. The 
sensitivity is maximum at m x ps 500 GeV, where cross 
sections as small as (av) ps 10~ 24 cm 3 s may be probed 
by IceCube. By observing the Milky Way halo, IceCube 
has already constrained the value of this annihilation 
cross section to be (av) < 10 -22 cm 3 s _1 for a DM mass 
of about 1 TeV [37]. We expect that the sensitivity 



obtained from observing the Virgo galaxy cluster will 
improve the above limit by about one order of magnitude 
if no detection is obtained. 

Han et al. [88], recently found evidence of extended 
gamma ray emission from the Virgo cluster, and the 
limit on the annihilation cross section that they obtained 
is (av) ps 10 _25 cm 3 s _1 . Although, in principle, this 
channel is observable at IceCube, we find that IceCube 
does not have the sensitivity needed to test this claim. 
Note that after the publication of the first version of this 
paper, it was found by several groups that the extended 
gamma ray emission in the Virgo cluster is due to the 
presence of new gamma-ray sources and not due to DM 
annihilation [120, 121]. 

The sensitivity that can be obtained by observing 
cascades is shown in Fig. 4 (top right panel). As 
KM3NeT is still under development, we show the 
constraints using two plausible choices for its angular 
resolution. We find that the sensitivity obtained from 
cascades is almost comparable with that obtained from 
muon tracks. The best sensitivity is achieved around 
a DM mass of around 20 GeV where a sensitivity to 
(av) « O(10 -25 cm 3 s _1 ) is reached, representing an 
order of magnitude improvement over the best sensitivity 
obtainable by observing tracks. 



B. XX vv 

In Fig. 4 (second from top and left panel), we show 
the expected sensitivity to self annihilation cross section 
(av) for XX ~~ * for m x in the range of (100 GeV — 
100 TeV), by detecting track-like events. The sensitivity 
is strongest at m x ps 500 GeV where the annihilation 
cross sections larger than (av) ps 0(lO _24 cm 3 s _:l ) can be 
probed. Due to the presence of substructures, we again 
get a sensitivity which is stronger by about one order 
of magnitude than the constraint obtained by IceCube 
when looking for this annihilation signal at the Milky 
Way Galactic halo [37]. Since the spectra of the signal 
and background are very different in this case, we expect 
that a much better sensitivity can be achieved due to 
a spectral analysis by the IceCube collaboration for the 
same exposure. 

We now consider cascade signals from this annihilation 
channel, in a DeepCore-like low energy extension in 
KM3NeT. For DM masses between 30 GeV and 10 TeV, 
the projected sensitivity is shown in Fig. 4 (second 
from top and right panel). In the low DM mass 
range, sensitivity to annihilation cross sections (av) ps 
O(10 -25 cm 3 s _1 ) can be reached. As can be seen from 
the plot, KM3NeT will have a unique opportunity to 
probe this part of the parameter space if it employs a 
low energy extension. 
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c. xx -> tt 

We now look at the sensitivity that can be obtained to 
the (av) for DM annihilation via \X ~^ ^t, by detecting 
track-like events. We show the sensitivity to (av), for m x 
in the range (1 TeV - 100 TeV), in Fig. 4 (third from top 
and left panel). Annihilation cross sections larger than 
(av) « 0(lO _23 cm 3 s _1 ) can be probed for DM mass in 
the range (1 TeV — 10 TeV). The constraints weaken for 
DM masses heavier than 10 TeV. 

Observation of cascades in a DeepCore-like low energy 
extension in KM3NeT give similar sensitivity in the 
(ITeV — 10 TeV) mass range. The expected sensitivity 
is shown in Fig. 4 (third from top and right panel). 
Annihilation cross sections (av) ~ O(10~ 23 cm 3 s _1 ) may 
be probed using the Virgo galaxy cluster. 

D. xx W vvvv 

In Fig. 4 (bottom left panel), we show the expected 
sensitivity to (av) for \X vVvV, for m x in the range 
of (100 GcV - 100 TeV), by detecting track-like events. 
The strongest sensitivity is achieved at m x sa 500 GeV, 
where the annihilation cross sections larger than (av) ss 
0(lO _25 cm 3 s~ 1 ) can be probed. Due to the presence of 
substructures, we again get a sensitivity which is stronger 
by about three orders of magnitude than the constraint 
obtained when assuming only an NFW profile. 

We now consider cascade signals from this annihilation 
channel, in a DeepCore-like low energy extension in 
KM3NeT. For DM masses between 10 GeV and 10 TeV, 
the projected sensitivity is shown in Fig. 4 (bottom 
right panel). In the low DM mass range, sensitivity to 
annihilation cross sections (av) « 0(1O -25 cm 3 s _1 ) can 
be reached. KM3NeT will have a unique opportunity 
to probe this part of the parameter space, which is not 
accessible by tracks, if it employs a low energy extension. 

Neutrino telescopes have not searched for neutrinos 
from the annihilation channel \X vvvv, but as we 
show in Fig. 4 (bottom panels), the constraints obtained 
in this channel can be quite promising. In particular [103] 
predicts enhanced emission in neutrinos, with (av) ~ 
10 _24 cm 3 s _1 in galaxies. The velocity dependence of 
the cross section in this model will reduce the cross 
section in galaxy clusters, but we believe that, besides 
the Milky Way and dwarf galaxies, galaxy clusters may 
also offer an important test for this model due to the 
strong substructure enhancement. 

E. Comparison with limits from other experiments 

In this section, we compare the sensitivity that 
IceCube and KM3NeT can obtain by observing clusters 
of galaxies with limits obtained from other experiments. 
The main annihilation products observable in an indirect 
DM detection experiment are neutrinos and gamma rays. 



Both of these signals are not deflected by galactic or 
extragalactic magnetic fields and come directly into the 
detector from the source. For the muon track signal, 
the IceCube results are for 276 days of live time and in 
the 22 string configuration. To compare our calculations 
with the same exposure in [37], we use 276 days as the 
time of observation in Eqs. (19) and (20) and take l/4th 
of the number of neutrino events to mimic the 22-string 
detector. For the cascade signal, we show our results for 
a 276 days of livetime of the detector and we take the 
volume available for cascades as 0.02 km 3 . 

We first show the various constraints on the \X 
/j + fi~ channel on the left panel in Fig. 5. This channel 
can be detected by both neutrinos and gamma rays. The 
decay of the muons produce electrons which can produce 
gamma rays via inverse Compton and bremsstrahlung. 
These energetic electrons can also produce synchrotron 
radiation in the galactic and extragalactic magnetic fields 
but the synchrotron radiation is at a lower frequency. 
The neutrinos are produced in this channel via the decay 
of the muon. 

At low DM masses (m x < 100 GeV), the constraints 
from Fermi-LAT using dwarf spheroidal galaxies are the 
most stringent [27]. This limit weakens for DM masses 
above a few hundred GeV. For higher DM masses, the 
limits on the DM annihilation cross section to muon 
pairs are obtained from the observation of the Milky 
Way Galactic halo by IceCube [37]. In the same figure, 
we show the improvements that can be obtained by 
observing galaxy clusters using neutrino telescopes. It is 
clear that galaxy clusters provide an order of magnitude 
more sensitivity compared to other sources. 

For this channel, we predict that the tracks observed 
in IceCube will give the best sensitivity for DM masses, 
m x > 1 TeV. KM3NeT, augmented with a DeepCore-like 
low energy extension, will provide increased sensitivity 
to this annihilation cross section at DM masses, m x < 
100 GeV. Although the expected sensitivity is weaker 
than the limits obtained from Fermi-LAT observation of 
dwarf galaxies, it will be an important complementary 
test, as the neutrino observations are less dependent on 
the central density profile. 

We show the various constraints on the \X ~ ^ v v 
channel on the right panel in Fig. 5. This indirect 
detection channel can only be detected by neutrino 
telescopes and it has no signatures in any other DM 
indirect detection experiment. 

For DM masses m x > 100 GeV, limits on the DM 
annihilation cross section to neutrino pairs are obtained 
from the observation of the Milky Way Galactic halo 
by IceCube [37]. In the same figure, we show the 
improvements that can be obtained by observing galaxy 
clusters using neutrino telescopes. Observation of the 
galaxy clusters by neutrino telescopes shall give an order 
of magnitude improvement over the existing constraints. 

For this channel, we predict that the observation of 
muon tracks in IceCube will give the best sensitivity 
above DM mass, m x w 300 GeV. KM3NeT augmented 
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with a DeepCore-like low energy extension will provide 
the best sensitivity to this annihilation cross section at 
DM masses m x < 1 TeV. For DM masses < 100 GeV, 
KM3NeT, with a low-energy extension, can reach 
annihilation cross sections of the order of 10 _24 cm 3 s _1 
while observing cascades in the detector. 



V. CONCLUSION 

In this paper, we have considered observation of 
galaxy clusters by neutrino telescopes and discussed 
the improvements that can be made over the existing 
limits. Recent high resolution computer simulations 
of galaxy clusters predict a large enhancement in the 
annihilation flux due to DM substructures. We take 
the substructure contribution into account and predict 
the neutrino flux from a typical galaxy cluster. We 
find that the sensitivity that can be obtained using 
galaxy clusters should improve the existing constraints 
by more than an order of magnitude. Our results should 
therefore encourage the IceCube collaboration to look at 
galaxy clusters, as an extension of their work on dwarf 
galaxies [28]. 

Due to the extended nature of the DM substructure 
profile (see Fig. 2), nearby galaxy clusters like Virgo 
should appear as extended sources at neutrino telescopes. 
We find that the optimal angular window around a galaxy 



cluster that maximizes the signal-to-noise ratio has a 
radius ps 2° (see Fig. 3). 

An order of magnitude improvement over the IceCube 
sensitivity is expected if KM3NeT deploys a low energy 
extension (like DeepCore in IceCube) in their telescope, 
which would allow for a full-sky observation with good 
pointing using cascades. This has the potential to open 
the (10 GeV — 100 GeV) DM mass range to neutrino 
astronomy, and improve existing constraints by an order 
of magnitude. We hope that these promising results will 
encourage the KM3NeT collaboration to investigate the 
possibility of deploying a low energy extension to their 
telescope and improve the reconstruction of cascades (see 
right panels in Fig. 4). 

We looked at the \X ~^ annihilation channel 

and predicted an order of magnitude improvement 
over the current constraints (see left panel in Fig. 5). 
Although this bound turns out to be weaker than 
the bound on the annihilation cross section given by 
Fermi-LAT while observing dwarf spheroidal galaxies, 
we emphasize that the large angular resolution of 
the neutrino telescopes make the result more model- 
independent than that obtained by Fermi-LAT. We have 
predicted that the improvement in sensitivity to the 
annihilation cross section in this channel will allow us to 
probe cross sections (av) > (10~ 24 — 10 _22 )cm 3 s _1 for 
DM masses in the range (1 GeV — 10 TeV) for 10 years 
of observation by a km 3 neutrino telescope. 
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We have also looked at the \X y v channel and 
predicted that the observation of galaxy clusters will 
constrain the annihilation cross section in this channel 
by an order of magnitude over the existing limit obtained 
by IceCube while observing the Milky Way Galactic 
halo (see right panel in Fig. 5). This annihilation 
channel is unique as it has no signal in any other 
DM indirect detection experiment. We predicted that 
the improvement in sensitivity to the annihilation cross 
section in this channel will allow us to probe cross 
sections (cry) > (1CT 24 — 10 _22 )cm 3 s _1 for DM masses 
in the range (1 GeV— 10 TeV) for 10 years of observation 
by a km 3 neutrino telescope. 

We considered the \X ~^ tt annihilation channel, 
which is expected to be very important for a heavy 
fermionic DM particle. We have predicted that the 
improvement in sensitivity to the annihilation cross 
section in this channel will allow us to probe cross 
sections (av) > 10 _22 cm 3 s _1 for 10 years of observation 
by a km 3 neutrino telescope. 

We finally considered the \X ~~ ^ vvvv channel and 
predict that the sensitivity that can be obtained using 
neutrino telescopes may be able to probe the annihilation 
cross sections required in models which aim to solve 
various small-scale problems in ACDM. 

Although we have performed our calculations for 
the Virgo galaxy cluster, we expect that neutrino 
telescope observation of a properly chosen galaxy cluster 
(after taking into consideration backgrounds and various 
detector systematics in more detail) will improve the 
limits on the annihilation cross section by an order of 
magnitude in almost all annihilation channels. We must 
emphasize that the biggest uncertainty in this result 



comes from the ~11 orders of magnitude extrapolation 
in the minimum DM substructure mass that is used to 
calculate the DM substructure profile. As a consequence 
of this extrapolation of the minimum substructure mass, 
the boost factor that can be obtained in a galaxy cluster 
due to the presence of substructures can vary by a factor 
of ~20. Unless simulations improve their resolution 
dramatically, this will remain an inherent assumption in 
any DM indirect detection experiment observing galaxy 
clusters. 

All things considered, we hope to have conveyed 
the usefulness of observing galaxy clusters at neutrino 
telescopes for studying DM. In particular, how good 
reconstruction of cascades can lead to significant 
improvements in sensitivity. We hope that the IceCube 
and the KM3NeT collaborations will consider our results 
and make the required improvements in their analyses 
and detectors to make this possible. 
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